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The nonrelativistic standard model for a continuous, one-parameter diffusion process in position 
space is the Wiener process. As well-known, the Gaussian transition probability density function 
(PDF) of this process is in conflict with special relativity, as it permits particles to propagate faster 
than the speed of light. A frequently considered alternative is provided by the telegraph equation, 
whose solutions avoid superluminal propagation speeds but suffer from singular (non-continuous) 
diffusion fronts on the light cone, which are unlikely to exist for massive particles. It is therefore 
advisable to explore other alternatives as well. In this paper, a generalized Wiener process is 
proposed that is continuous, avoids superluminal propagation, and reduces to the standard Wiener 
process in the non-relativistic limit. The corresponding relativistic diffusion propagator is obtained 
directly from the nonrelativistic Wiener propagator, by rewriting the latter in terms of an integral 
over actions. The resulting relativistic process is non-Markovian, in accordance with the known fact 
that nontrivial continuous, relativistic Markov processes in position space cannot exist. Hence, the 
proposed process defines a consistent relativistic diffusion model for massive particles and provides 
a viable alternative to the solutions of the telegraph equation. 
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I. INTRODUCTION 



It is well-known for a longtime that the nonrelativistic 
diffusion equation (DE) [l|,H HI, [| 
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is in conflict with the postulates of special relativity. A 
simple way to see this is to consider the propagator of 
Eq. |T]), which for d = 1 space dimensions, is given by 
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[Heaviside's theta-function is defined by <d(t) = l,t > 
and 0(i) = O, i<0]. This propagator @ repre- 
sents the solution of Eq. (JTJ) for the initial condition 
p(0,x) — S(x — xo). That is, if X{t) denotes the path 
of a particle with fixed initial position X(0) = xq, then 
pit, x\xq)c\x gives the probability that the particle is 
found in the infinitesimal volume clement [x,x + dx] at 
time t > 0. As evident from Eq. for each t > there 
is a small, but non- vanishing probability that the particle 
may be found at distances \x — xq\ > ct, where c is the 
speed of light. 

This problem has attracted considerable interest dur- 
ing the past 100 years (see, e 
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and references therein). Nonetheless, it seems fair to say 
that a commonly accepted solution is still outstanding. 
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Apart from the profound theoretical challenge of develop- 
ing a consistent relativistic diffusion theory, several prac- 
tical applications exist, e.g., the analysis of data from 
high energy collision experiments [IE EE HE HE HI] ■ In 
this context, an often considered alternative to Eq. ( fT|) is 
given by the t eleg raph equation (TE) [E EE EE Il4 Il5l 

EE SI ISIS [II 
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where Td is an additional relaxation time parameter. As 
e.g. discussed by Masoliver and Weiss [101 or Abdel- 
Aziz and Gavin [IE], the diffusion fronts described by 
Eq. ([3]) spread at finite absolute speed a = ±(D /t^) 1 / 2 ; 
cf. Fig. 1 of Ref. [HI]. In particular, by fixing Cd = c, 
the solution of Eq. ([3]) defines a simple relativistic, non- 
Markovian [EH diffusion process without superluminal 
propagation. The 'nonrelativistic limit' corresponds to 
letting Cd — > oo in Eq. ([3]), which gives Td — *■ and thus 
leads back to Eq. (P). 

However, apart from these appealing properties the so- 
lutions of Eq. (|3|) exhibit a peculiar feature that is rather 
unexpected from a simple diffusion theory: The diffu- 
sion fronts are given by running <5-peaks whose inten- 
sity decreases as time grows; cf. Eq. (32) in Ref. [To| . 
The appearance of such singularities can be understood 
by the fact that the one-dimensional TE ([3]) can be 
derived from a so-called persistent random walk (RW) 
model [E HE HE HE HE]- This model assumes that a 
particle moves with constant absolute velocity between 
neighboring lattice points. At each lattice point, the 
particle is either back-scattered or transmitted, with 
the transmission probability being larger than the back- 
scattering probability (persistence). As Goldstein § 
showed for the one-dimensional case, the continuum limit 
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of this model leads to the TE (see Bogufia et al. [H| for 
a generalization to higher space dimensions). 

Hence, in contrast to the ordinary DE (JXJ) , the 
TE ([3]) relies on asymmetric transition probabilities, 
which among others lead to a non-vanishing probability 
concentration at the diffusion fronts. From the practi- 
cal point of view, the concept of persistent diffusion can 
certainly be useful in some situations, e.g., if one aims 
at descri bing the propagation of photons in thin slabs 
or foams [ill l36l 13?! \3&[ . In general, however, it seems 
unlikely that the relativistic generalization of the sim- 
ple nonrelativistic diffusion theory would require concep- 
tual modifications as severe as the introduction of per- 
sistence. In particular, for massive particles <5-peaked 
diffusion fronts are unlikely to occur in nature, because 
such fronts would imply that a finite fraction of parti- 
cles carries a huge amount of kinetic energy (much larger 
than mc 2 ). One may, therefore, wonder if it is possible to 
construct alternative, non-persistent diffusion processes, 
which are consistent with the basic requirements of spe- 
cial relativity, and converge to the solutions of the non- 
relativistic DE ((T|) in a suitable limit case. This is the 
question the present paper aims to deal with. 

Before outlining our approach, a general remark might 
be in order. Usually, a diffusion theory intends to pro- 
vide a simplified phenomenological description for the 
complex stochastic motion of a particle in a background 
medium (e.g., on a substrate 5, 30, 33ll34l.l35l| or in a heat 
bath [20|V Thus, there exists a preferred frame, corre- 
sponding to the rest frame of the substrate (or, more gen- 
erally, the center-of-mass frame of the interaction sources 
causing the stochastic motion). It is therefore not expe- 
dient to look for Lorentz or Poincare invariant spatial 
diffusion processes (cf. Sec. 5 of Montesinos and Rov- 
elli [39l | ) . Accordingly, we focus here on discussing simple 
diffusion models that comply with the basic requirements 
of special relativity in the rest frame of the substrate. 

Very often, diffusion models are based on time evolu- 
tion equations like the DE ([T]), derived either from an 
underlying microscopic model or as an approximation to 
a more complicated phenomenological model [l], H, 0, H3, 
miH. Unfortunately, apart from the TE ©, the con- 
struction of relativistic DEs turns out to be very cumber- 
some. Among other things, this can be attributed to the 
fact that standard limitingjsrocedures, as those leading 
from discrete RW models 2] to Eq. ([I}, typically result 
in superluminal propagation (cf. the discussion in Sec. HT1 
below). Therefore, a different approach is pursued in the 
present paper: Instead of aiming at a relativistic evolu- 
tion equation in the first place, we try to obtain a con- 
tinuous, relativistic generalization of Eq. $2$ by focussing 
on the structure of diffusion propagators. 

The paper is organized as follows. First, some ba- 
sic properties of Markovian diffusion models for massive 
particles are reviewed (Sec. |TT|. In particular, it will 
be demonstrated in detail why it is impossible to con- 
struct a continuous relativistic Markov model in position 
space (Sec. Ill Cj) . Subsequently, we propose a simple non- 



Markovian generalization of the nonrelativistic Wiener 
process (Sec. IIIip . The specific functional shape of the 
relativistic propagator arises naturally, if one rewrites the 
Wiener propagator ^ in terms of an action integral. The 
paper concludes with a summary of the main results. 

II. MARKOVIAN DIFFUSION MODELS 

In this section, we consider a standard Markovian RW 
model that provides the basis for the subsequent consid- 
erations. First, the underlying equations and assump- 
tions are summarized. Then it is discussed why the dis- 
crete RW model may lead to a useful nonrelativistic con- 
tinuum model, but fails to do so in the relativistic case. 
For clarity, we focus on the one-dimensional case d = 1 
in this part [57| . 

A. Basic equations and assumptions 

Consider a massive particle located at some known po- 
sition X(0) — xq at the initial time to = 0. Assuming 
that the particle performs random motions, we are inter- 
ested in the distribution of the particle position X(t) = x 
at time t > 0. To obtain a simple RW model for the 
stochastic dynamics of the particle, we start from the 
formal identity 

X(t) = x + / ds V{a), (4) 
Jo 

where V(s) is the velocity. Next we assume that the 
motion of the particle can be viewed as a sequence of free 
ballistic flights, interrupted by (N — I) scattering events. 
Mathematically, this idea is formalized by introducing 
the discretization T/v = {to, . . . , ijv} of the interval [0, t], 
such that 

= tg < t\ < . . . tjV — 1 < t y\r = t. 

Defining the times-of-free-flight by 

n := t t - ti-i, i = l,...N, 
we may discretize Eq. ((4]) by using the approximation 

JV 

X(t) = x + '^2nv i . (5a) 

i=l 

Here, Uj = V(U) is the constant velocity of the particle 
between tj_i and ti, i.e., before the ith collision. 

Within the standard Markovian approach, it is com- 
monly assumed that the velocities (vi,V2, ■ ■ ■) can be 
viewed as independently distributed random variables 
with a given PDF 

JV 

fN(vi,-.-,v N ;{Ti}) = Y[f(vi;Ti). (5b) 
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For any fixed partition {Tj}^^...^, this uniquely deter- 
mines the transition PDF for the position coordinate as 



p N (t, x\x ) 



N 

n j dvif^Ti 



S(x~X(t)), (5c) 



where X(t) is the discretized path from Eq. ((Saj). 
Throughout, we abbreviate PNit, x\xq) = PN(t, x\to, xq) 
if the initial time is fixed as to = 0. 

The continuum limit is recovered by letting N — > oo, 
such that r, — > and 



N 



}^Tj=t-t =t. 



(6a) 



The transition PDF in the continuum limit can then be 
expressed as 



p(t, x\x ) = lim p N (t,x\x ) 



(6b) 



Yet, the continuous Markov model ([6]) is useful only if 
this continuum limit yields a nontrivial result. As a fa- 
miliar example, we will recall how the (marginal) velocity 
distribution /(v;Tj) is to be chosen in order to recover 
the nonrelativistic diffusion propagator @ from Eqs. ([5]) 
and (J6j) . Thereby, it will become evident why this Marko- 
vian approach fails to produce a useful continuum model 
in the rclativistic case. 



B. Nonrelativistic diffusion model 



In the nonrelativistic case, we choose 



/(Wi;r ^(i^) 1/2exp (; 
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This can be interpreted as a Maxwellian velocity distri- 
bution with an inverse (effective) temperature 



a = (fc B T,r 



2mD 



(8) 



where m denotes the particle rest mass, and fee the Boltz- 
mann constant. Inserting Eq. ([7]) into Eq. (|6b[) gives 



p(t, x\xq) = lim 

N^oc 
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FIG. 1: Spreading of the Gaussian PDF p(t,x) — p(t,x\0) 
from Eq. ((2} at different times t, where t is measured in units 
of D/c 2 . At initial time t — 0, the PDF corresponds to a 
5-function centered at the origin. 



it is straightforward to show that the rhs. of Eq. © 
is equal to the nonrelativistic diffusion propagator from 
Eq. ©, and, hence, defines a Wiener process [43| . Then, 
fixing the element of the functional integration by 



V[V] 



lim 

Af^oo 



1/2 



(10) 



one finds the following functional integral representation 
of the nonrelativistic diffusion propagator 



P{t, x\x ) 



V[V] exp 

S[ (x - X ) 



f dsV(s] 
Jo 
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Very often, path integral representations of diffusion 
propagators [40l . |44| are formulated in terms of curves 
X(t) in position space, whereas Eq. (fTTj) expresses the 
diffusion propagator ^ by means of a functional inte- 
gral over the space of the velocity curves. This is more 
convenient with regard to the subsequent discussion. 

The nonrelativistic PDF (|2"TTTj) is depicted in Fig. [TJ 
For later use, we mention that the mean square displace- 
ment of the Wiener process is given by [l| 

([X(t)-x ] 2 ) := J dx (x-x ) 2 p(t,x\x ) = 2Dt. 

With regard to the subsequent discussion, however, it 
is important to notice that the continuum limit Tj — > 
implies that the effective temperature Tj of the Maxwell 
distribution ([7]) diverges, Tj — * oo; i.e., the second mo- 
ment (v 2 ) of the velocity distribution /(fjjTj) tends to 
infinity in the continuum limit. As discussed in the next 
part, this divergence is essential for ensuring that those 
scenarios, where a particle leaves its initial position, re- 
tain a non- vanishing probability in the continuum limit. 
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C. Relativistic case 

The superluminal propagation of particles can be 
avoided in the discrete model by considering appropriate 
velocity distributions Tj) with finite support [— c, c]. 
To keep subsequent formulae as simple as possible, we 
shall use natural units with c = 1 (even though, c will 
be reinstated occasionally in figure captions or to clarify 
units). Moreover, for simplicity, we focus on the case of 
an equally spaced partition 



n = r = t/N, Vi = l,...,N, 

and write /(i>i) instead of t). 



(12) 



1. Time-discrete Markov model 

The relativistic generalization of the one-dimensional, 
nonrelativistic Maxwell-distribution is given by the 
one-dimensional modified Juttncr function 20, 45] 



f(Vi) = 



7fa) 2 
2^o(x) 



exp[-xi(vi)] 9(1 -|«,|), (13) 



where x is a (dimensionless) temperature parameter, and 
-Ko(x) a modified Bessel-function of the second kind [46| . 
Furthermore, 

E = nvy(v), p — mrf{v) v, j(v) — (1 — w 2 ) -1 ^ 2 

denote the relativistic kinetic energy, momentum and 7- 
factor, respectively. Analogous to the Maxwell distribu- 
tion in the nonrelativistic case, the PDF (fl~3"|) is distin- 
guished by the fact that it is conserved in relativistic 
elastic binary collisions of classical particles [2(| ■ 

From the nonrelativistic limit, i.e., by expanding the 
exponent of Eq. (| 13|) for u? -C 1 and comparing with 
Eq. ((TJ, one can identify 

rn t 

X = m(3 



knT 2D 



Upon inserting the PDF (fT3|) into Eq. (|5c|) . one obtains 
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In contrast to the nonrelativistic case, it is very difficult 
or, perhaps, even impossible to evaluate the integral (|14|) 
analytically for N > 1. The propagator for the trivial 



case N = 1 is given by 
-1 



pi(t,x\x ) 
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Before discussing multiple scatterings N > 1 we note 
that, in principle, there are two ways of reading the 
Eq. (fT2|) : If one fixes t and increases N, then this corre- 
sponds to taking the continuum limit r — > 0. Alterna- 
tively, one can consider the time-of-free-flight r as being 
fixed; then the discretization (|12p just means that t is 
given by multiples of r. In the remainder of this sub- 
section, we shall adopt the latter point of view, because 
this leads to a nontrivial time-discrete relativistic Markov 
model. By contrast, as discussed afterwards, the contin- 
uum limit t — ► fails to produce a useful diffusion model 
in the relativistic case. 

However, assuming for the moment that r is fixed, the 
probability density at time t = Nt for a particle with 
initial condition xq = is given by 

p N {x) :=pn(Nt,x\0). 

Numerical results for different values of and two dif- 
ferent parameter values x — rnc 2 /(k^T) or D = r/(2x), 
respectively, are depicted in Fig. [5J As evident from the 
diagrams, superluminal particle propagation does not oc- 
cur anymore, due to the choice of the velocity PDF (|13jl . 
Moreover, as the comparison with Fig. [T] shows, at large 
times t ^S> t the results start to look similar to those of 
the corresponding nonrelativistic diffusion process. 

Furthermore, for a strongly relativistic velocity PDF 
with X i$ 1) see Figs. [2] (b), we also observe peculiar peaks 
due to the discrete time structure of our model. Some- 
what miraculously, such effects are absent in the nonrel- 
ativistic model (cf. Sec. Ill B[) . but this is attributable 
to the fact that the Maxwellian (Gaussian) PDF repre- 
sents a stable distribution JS?]]; i.e., it is form-invariant 
under consecutive foldings as those in Eq. (f5c]l . In gen- 
eral, however, this type of invariance does not hold any- 
more if the velocity distribution is non-Gaussian and of 
finite support. In the latter case, characteristic relax- 
ation patterns (peaks, dints, etc.) may arise which reflect 
the space/time discretization scheme underlying the RW 
model. 

For instance, the appearance and positions of the peaks 
in Figs. [2] (b) can be readily understood by considering 
the ultra-relativistic case x — ^ 0, where the velocity dis- 
tribution (fl"3|) converges to two 5-peaks at ±c: Let us 
assume that an ultra-relativistic particle starts at time 
t = at position xq = 0. Then at time t — r we have 
p(x, t|0, 0) = [S(x + t) + 5(x — r)]/2, corresponding to 
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FIG. 2: Transition PDFs p N (x) = p N (Nr, x\0, 0) for the rel- 
ativistic time-discrete RW model from Sec. Ill C II The PDFs 
were numerically calculated at five different times t = Nt for 
two different parameter values \ = rnc 2 /(k^T) — r/(2D). 

(a) Weakly relativistic velocity distribution with x > 1- 
The numerical solutions look qualitatively similar to a (non- 
relativistic) Gaussian, but have a finite support [— ct, ct]. 

(b) Strongly relativistic velocity distribution with x < 1- The 
peaks reflect the discrete time steps of the model, and are as- 
sociated with the backward/forward-scattering of the fastest 
particles; cf. explanation in the text. 



two peaks of equal height. Analogously, at t — 2r one 
finds 2t|0, 0) = [5(x + 2r) + 2S(x) + 5(x - 2r)]/4 
corresponding to three peaks, etc.. 

From the general point of view, however, it would be 
desirable to also have continuous, analytically tractable 
diffusion models at one's disposal which do not exhibit 
such peculiar relaxation effects and discretization signa- 
tures. Such a model will be proposed in Sec. IIIH Before, 
it is worthwhile to study in more detail why the contin- 
uum limit t — ► does not provide a useful Markov model 
in the relativistic case. 



tic analogon of the path integral (flTj) as 
p(t,a;|a;o) = / V[V] exp 



2D 



S[ (x-xo)- J o dsV(s) ], (16a) 



where 



2K, 



7(^) 2 1 
o[t/(2DN)} j 



(16b) 



However, as we shall see immediately, this implies the 
trivial result 



p(t, x\xq) — 5(x — xq), V t > 0. 



(17) 



The proof can be performed most easily by means of the 
Central Limit Theorem (CLT) [H, H§| . For this purpose, 
consider the mean velocity along a path, defined by 



X(t) -logglr^ GUl-r^ 

V N {t):= = j}^ T ^ = n!^ Vi - 

i=i »=i 



The random variables Vi are identically, independently 
distributed, with mean (vi) = and variance (vf \ < oo. 
For this case, the CLT asserts that the distribution of 
the random variable Zjv(t) := VNVwit) converges to a 
Gaussian. Hence, in the continuum limit N — > oo, the 
mean velocity Vpj(t) goes to zero with probability one, 
which means that the particle effectively remains at its 
initial position, i.e., p(t, x\xq) = S(x — xq) as anticipated 
above. 

One may therefore conclude that a nontrivial continu- 
ous, relativistic Markov process in position space cannot 
exist [58| . The main reason is that for any relativistic 
velocity distribution the variance is bounded (yf) < oo 
because the values Vi are bounded, and thus the CLT 
applies. This also means that, in principle, it is not very 
promising to consider relativistic path integrals of the 
Markovian type (|16ap . In particular, in view of the close 
analogy between diffusion and quantum propagators in 
the nonrelativistic case, this seems to leave very little, 
if any, room for formulating a relativistic particle quan- 
tum mechanics in terms of imaginary time Markov path 
integrals 44 1. 



III. NON-MARKOVIAN CONTINUOUS 
DIFFUSION MODELS 



2. Shortfall of the relativistic continuum limit 

Keeping t fixed and taking the limit N — > oo of the 
time-discrete propagator (|14[) . one can write the relativis- 



The preceding discussion has shown that one has to 
abandon the Markov property if one wishes to define a 
continuous diffusion model in position space that avoids 
superluminal particle propagation speeds. In this section, 
we identify a simple non-Markovian generalization of the 
Wiener propagator Q that is confined to the light cone 
and exhibits continuous diffusion fronts (in contrast to 
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the solution of the telegraph equation). As we shall see, 
the specific functional shape of the proposed relativistic 
propagator naturally arises from the nonrelativistic case, 
provided one rewrites the transition PDF of the Wiener 
process in terms of an integral over actions. 

We generalize the discussion to an arbitrary number 
of space dimensions d from now on. Abbreviating x — 
(t,x) and xo — (to,Xo), the transition PDF of the d- 
dimensional Wiener process is given by 



p(x\x ) = N d exp 



(x - x ) 2 
AD{t-t ) 



(18) 



where the normalization constant Nd is a function of the 
time difference (ii — to) > 0. Analogous to the one- 
dimensional case from Sec. Ill B[ the Wiener propaga- 
tor (Tl8|) could be expressed via a standard Markovian 
path-integral construction. But again, this type of repre- 
sentation is of no use with regard to the relativistic gener- 
alization for the same reasons as discussed in Sec. Ill C 21 
Hence, we shall next look for another formal representa- 
tion of Eq. (fl8|) , that can be consistently transferred to 
the relativistic case. 

For this purpose, consider a particle traveling from the 
event xq — (to,Xo) to x — (t,x). Assume that the par- 
ticle can experience multiple scatterings on its way, and 
that the velocity is approximately constant between two 
successive scattering events. Then the total action (per 
mass) required along the path is given by 



dt' v(t') 2 



(19) 



where the velocity v{t') is a piecewise constant func- 
tion. Clearly, the action becomes minimal for the de- 
terministic (direct) path, i.e., if the particle does not col- 
lide at all. In this case it moves with constant velocity 
v(t') = {x — xo)/(t — to) for all t' £ [to,t], yielding the 
minimum action value 



r — — \ 1 (X- Xq) 2 



(20) 



On the other hand, to match the boundary conditions 
it is merely required that the mean velocity equals 
(x — Xo)/(t — to). Consequently, in the nonrelativis- 
tic case, the absolute velocity of a particle may be- 
come arbitrarily large in some intermediate time interval 
[t',t"} C [to,t]. Hence, the largest possible action value 
is a + (x,x~o) = +oo. These considerations put us in the 
position to rewrite the Wiener propagator (TT8")) as an in- 
tegral over actions: 



p(x\x Q ) = N' d 



da exp 



v 2£>)' 



(21) 



The important advantage of this representation lies in the 
fact that it may be generalized to the relativistic case in a 
straightforward manner: One merely needs to insert the 
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FIG. 3: Transition PDF p(t,x) = p(t,x\0) for the one- 
dimensional (d = 1) relativistic diffusion process (|23[) at dif- 
ferent times t (measured in units of D/c 2 ). 



corresponding relativistic expressions into the boundaries 
of the integral. 

The relativistic generalization of Eq. (fT9|) reads 



J t 



(22) 



Analogous to the nonrelativistic case, the relativistic ac- 
tion (|22p assumes its minimum a_ for the deterministic 
(direct) path from Xq to x, characterized by a constant 
velocity v(t') = (x — xq) /(t — to). One explicitly obtains 



a-{x, x ) = -[(t- t ) 2 -{x- xof 



I 1/2 



(23a) 



i.e., a_ is the negative Minkowski distance of the two 
space-time events x~o & n d x. The maximum action value 
is realized for particles moving at light speed, yielding 
a + = 0. Hence, the transition PDF for the relativistic 
generalization of the Wiener process reads 



p(x|a;o) = JV d \ exp 



a-(x,x ) 



2D 



(23b) 



if (x — Xq) 2 < (t — to) 2 , and p(x\xo) = otherwise, 
with a- determined by Eq. (|23ap . Explicit results for 
the normalization constant Nd with d = 1, 2, 3 are given 
in App. El 

Figure [3] shows the PDF p(t,x) = p(t,af|0,0) of the 
diffusion process (|23"]) for the one-dimensional case d = 1 
at different times t. The corresponding mean square dis- 
placement is plotted in Fig. 2] (dashed curve). 

The diffusion process defined by Eqs. (|2"3"]) is non- 
Markovian. This can be seen, e.g., by testing the nec- 
essary Chapman-Kolmogorov criterion 

p(t, x\t ,x ) = / d d x 1 p(t, x\t ± ,xi) p(ti,xi\x ). (24) 



If Eq. becomes invalid for some t\ € (to,t), then the 
process is non-Markovian. As one my easily check, e.g., 
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FIG. 4: Comparison of the mean square displacements 
(X 2 {t)), divided by 2Dt, for the one-dimensional (d — 1) non- 
relativistic Wiener process @ and its relativistic generaliza- 
tion from Eq. (|23[) with initial condition (to,xo) = (0,0). 



by numerically evaluating (|24p for some sample values 
(t,x;to,Xo) and some intermediate time t±, Eq. is 
indeed violated for the transition PDF (|23| . 

It is interesting to note that the PDF (|2"Tj) is a special 
case of a larger class of diffusion processes, defined by 



/•a+(x|x ) 

p w (x\x ) = Af w / daw(a), 

J a_ (x|xo) 



(25) 



where w(a) > is a weighting function, and Af w the 
time-dependent normalization constant. Equation (|21[) is 
recovered by choosing u>(a) = exp[— a/ (2D)]. It is, how- 
ever, worth mentioning that a very large class of func- 
tions w(a) yields an asymptotic growth of the spatial 
mean square displacement that is proportional to i, cor- 
responding to 'ordinary' diffusion. Moreover, Eq. ([23]) 
can also be used to describe superdiffusion or subdiffu- 
sion processes [5l| HH, [53| , whose asymptotic mean square 
displacements grow as i Q , a ^ 1 [53 ]. 

In particular, Eq. (|25p may be viewed as an 'un- 
conventional' non-Markovian path integral definition in 
the following sense: Physically permissible paths from 
xq to x have action values (per mass) a in the range 
[a_, a+]. Grouping the different paths together according 
to their action values, one may assign to each such class of 
paths, denoted by C(a; x, xq), the statistical weight w(a). 
The integral ([25]) can then be read as an integral over 
the equivalence classes C(a;x,Xo) and their respective 
weights w(a). The nonrelativistic Wiener process cor- 
responds to the particular choice w(a) = exp[— a/(2D)]; 
hence, it is natural to define the relativistic generalization 
by using the same weighting function. 



IV. SUMMARY 



The search for a consistent description of diffusion 
processes in the framework of special relativity repre- 
sents a longstanding issue in statistical physics. Sev- 
eral different relativistic generalizations of the nonrela- 
tivistic diffusion equatio n fl] ) have been proposed in the 
past @, EE El ES El El ES El- The most promi- 
nent example is the telegraph equation ([3]). This sec- 
ond order partial differential equation (PDE) describes 
a non-Markovian diffusi on p rocess without superluminal 
particle propagation [5j, |30|, |3l|, |32j. However, the solu- 
tions of the telegraph equation exhibit singular diffusion 
fronts propagating at light speed c, and may therefore be 
less appropriate for the description of massive particles. 

To avoid such singularities, we pursue a complemen- 
tary approach in this paper; i.e., instead of trying to find 
a relativistic diffusion PDE we focus on the transition 
probability density function (PDF) of a relativistically 
acceptable diffusion model. As the main result, an an- 
alytic expression for the transition PDF is proposed in 
Eq. ([2"5|) . The simple relativistic diffusion process de- 
fined by this PDF avoids superluminal propagation, and 
reduces to the Wiener process ([2]) in the nonrelativis- 
tic domain. This analytically tractable relativistic diffu- 
sion model can be useful for practical calculations in the 
future (e.g., in high energy physics [IH, EB| and astro- 
physics [541]). 

Finally, it is worth emphasizing that, analogous to the 
solutions of TE but in contrast to the Wiener process, 
the proposed diffusion model is non-Markovian. This is 
not a drawback, as it is generally impossible to construct 
nontrivial continuous relativistic Markov processes in po- 
sition space (cf. discussion in Sec. Ill Cl and Ref. @). In 
view of the close analogy between diffusion and quantum 
propagators in the nonrelativistic case, it would be quite 
interesting to learn, if this also means that it is impossible 
to formulate a relativistic particle quantum mechanics on 
the basis of standard path integrals 0, HH . 
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APPENDIX A: CALCULATION OF THE 
NORMALIZATION CONSTANT N d 

We wish to express the normalization constant Md from 
Eq. (|23bp in terms of modified Bessel functions of the 
first kind and modified Struve functions [461 ] . Introducing 
z := x — Xq and u :— t — t$, we have to calculate 



Md = \ d d zG(u-\z\)\exp 



(u 2 -z 2 ) 



2)1/2 



2D 

Using spherical coordinates, we can rewrite this as 

m-i L„J(^ 2 -N 2 ) 1/2 



1 



Af d = O d d\z\ \z 



exp 



2D 



- 1 



(Al) 

where O d = 2-K d l 2 /T(d/2) is the surface area of the d- 
dimensional unit-sphere. It is convenient to split the in- 
tegral (|Alj) in the form 



(A2) 
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where 

M = O d / d|z| \z\ d - 1 



exp 



(u 2 - \z\ 2 ) 



2U/2 



2L> 



(A3) 



Next we substitute |z|= it sin 0, where cf> 6 [0,7r/2]. 
Then Eq. (|A"3|l takes the form 



/•t/2 

A/^ = u d O d ■ 
Jo 



cos sin 



d-l 



exp 



u cos < 
2D 



Based on this integral representation, M' d can be ex- 
pressed in terms of modified Bessel functions of the first 



kind I n , and modified Struve functions Lk [461 ] ■ and one 
finds 



utt [7 1 (x)+i-i(x)] 



2tt 



X 

,2 



[1 + (X - 1) exp(x)] 



(A4a) 
(A4b) 



u 3 ix IHx) + Lo(x)} - 2Li(x)} , (A4c) 



X 



where x = u/(2D). 



